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Abstract 

Effect of site-selective substitution of Ti in Cu on the thermal stability 
of CuTi alloy is investigated using classical molecular dynamics simula- 
tions with Embedded Atom Method potentials. It has been observed 
experimentally that melting point of all the naturally occurring stable 
phases of CuTi alloys do not show a definite trend with gradual increase 
in Ti concentration. To understand the phenomenon, super cells of CuTi 
alloy are constructed where Cu atom is substituted by Ti randomly and at 
selective sites. For random substitution, the melting point decreases lin- 
early with increase in Ti concentration. A non-monotonous dependence is 
seen when Cu atoms at selective sites are replaced by Ti. For a particular 
doping concentration, the melting point shows a wide range of variation 
depending on the order of atomic arrangement, and can be fine tuned 
by selecting the sites for substitution. The variations in melting points 
in different cases are explained in terms of the peak height, width and 
position of the corresponding radial distribution functions. Finally, it is 
verified that initial structures of the naturally occurring CuTi alloys are 
responsible for the non-definite trend in their melting points. 
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1 Introduction 



Identification of doping strategies is important to optimize the physical proper- 
ties of alloys . Usually dopant concentration is varied to tailor the properties. 
For a fixed concentration, position of the dopant atom in the host matrix is also 
a key factor as it determines the interaction between host and dopant atoms. 
Thus doping at the selective sites is a viable tool to achieve finer control over 
the properties of the alloy. Site selective doping has been realized experimen- 
tally and is also used to make an educated guess of the material properties. For 
example, spatially controlled doping of a set of impurity atoms (cobalt, silver, ni- 
trogen, and boron) in carbon nanotube has been achieved experimentally by ion 
implantation technique pQ. Lu2Si05 has been doped by Eu 3+ at preferential site 
by sol-gel chemistry method [2] . Site selective substitution of boron in carbon 
nanoribbons has been investigated in the framework of first-principles density 
functional theory. Electronic properties are obtained for different boron— boron 
arrangements and concentrations 3J. While majority of these reports are fo- 
cused on the electronic properties, role of site selective doping on the thermal 
stability of the alloy is rarely investigated. 

In the present work we investigate the role of site selective substitution of Ti 
in Cu on the melting point of CuTi alloy. CuTi is technologically significant due 
to its high strength and conductivity. It has high yield strength when hardened 
via precipitation, is widely used for direct brazing and joining of ceramics j4] 
and considered as a promising material for thermoelectric power generation [5] . 
There are three main aspects of studying variation in melting point of CuTi 
alloys. 

The cell volume of some alloys increases compared to that in its pure phase 
due to the larger radius of the dopant atoms. Melting point of these alloys 
decrease if the melt occupies larger volume than its solid phase. CuTi alloy 
belongs to this category and the melting point of CuTi alloy should decrease as 
Ti concentration increases. But the phase diagram of CuTi alloy shows decrease 
in melting point, from 1356 K to 1158 K, on alloying by Ti only below a certain 
concentration (20 % in Cu4Ti) [6]. Beyond this limit, stable phases are found at 
much elevated temperatures. Melting point gradually increases from 1158 K to 
1273 K as the Ti concentration increases for stable phases like Cu2Ti, Cu3Ti2, 
Cu4Ti3 and CuTi (50 %). Therefore linear dependence of melting point on 
Ti concentration is not observed. This anomaly in the melting curve can be 
explained by performing simulations on thermal stability of CuTi alloy having 
various arrangements of dopant atoms. 

Crystals of CuTi alloys can be prepared over a wide range of Ti concentration 
by heat treatment and quenching of the melt of pure Cu and Ti. Cu-Ti phases 
having different thermal stabilities can be obtained by changing the temperature 
and duration of heat treatment [7], [8]. Also, hardening of CuTi alloys with low 
Ti concentration, via aging is an important technique to enhance its ductility. 
Both the methods of preparation and processing, like age hardening of different 
phases of CuTi alloy, require knowledge of melting point in advance. 

There are issues related to precipitation during age hardening of alloys to 
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create high yield strength. In this process clusters of dopant atoms are created 
by cooling the heat treated alloy for forming mobile defects. Since these kinds 
of defects have large impact on the melting point, their creation requires prior 
information on the thermal stability of the alloy. 

The organization of the paper is as follows: In section^ we discuss classi- 
cal molecular dynamics method for studying the CuTi system. The one phase 
method for determining the melting point, the embedded atom method (EAM) 
potential and the different initial structures obtained through doping are intro- 
duced. In the section[3]we discuss the impact of random substitution, impurity 
atom cluster/micro-structure formation and site-selective substitution on the 
melting point of the alloy. The anomaly in the naturally existing stable phases 
of CuTi alloys is also considered. The variations in melting point are correlated 
with the changes in peak position, height and width of the radial distribution 
function (RDF) . Finally our results are summarized in section 01 

2 Computational Details 

2.1 Classical Molecular Dynamics simulation 

The melting points of the CuTi alloy are determined using classical molecular 
dynamics simulations [3] . We employ the parallel molecular dynamics simulation 
package DL_POLY [TU]. DL_POLY has been used earlier to model the melting 
process of the inter- metallic compound N13AI [TT]. We consider cubic super 
cells consisting of 11 x 11 x 11 conventional FCC unit cells, which corresponds 
to 5324 atoms. Three dimensional parallelopiped periodic boundary conditions 
are applied on the super cells to eliminate the surface effects and reproduce the 
bulk properties. The Berendsen isothermal-isobaric (NPT) ensemble is used to 
achieve constant temperature and pressure conditions [12] . The relaxation times 
for the thermostat and barostat are 1.0 and 3.0 ps, respectively, and the pressure 
is fixed at atm. We use the Verlet leapfrog scheme for integrating the Newton's 
equations of motion [13 with a time step of 0.001 ps for all cases. Simulations 
are done for a total of 10,000 time steps, where the first 4000 steps are used 
for equilibration and the remaining 6000 for statistical averaging. In order to 
ascertain that only one image of a particle interacts with another particle, the 
cut off distances for force calculations is chosen to be smaller than half the size 
of the supercell. A cutoff distance of 10A is used for all the simulations. 

2.2 Determination of melting point using one phase method 

The melting point of the CuTi systems is determined using the one-phase 
method [14] wherein the melting point is identified from sharp increase in atomic 
volume, diffusion coefficient and energy as temperature is varied. Initially the 
temperature is increased by 100K until the alloy melts and the range within 
which the jumps occur is observed. The structure obtained after equilibration 
at a certain temperature is used as the starting configuration for the subsequent 
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Figure 1: Radial distribution function (RDF) of Cu before and after melting. 

higher temperature. The melting point for the specific alloy is next obtained 
by increasing the temperatures by 10K and observing the jumps in the above 
mentioned quantities. Also, as the long range order is lost once melting occurs, 
the height of the first peak of RDF decreases drastically and RDF m 1 after 2-3 
small humps. The RDF describes how the atomic density varies as a function 
of distance from the atom and is given by g(r) — n(r) / (Airr 2 Ar p) where n(r) is 
the number of atoms at a distance r within a shell of thickness Ar and p is the 
average number density. The heights of the peaks signify the number of first, 
second, etc., nearest neighbours. Similarly the position of RDF peaks reflects 
the neighbour distances. Earlier, a similar study of alloying bcc Fe with Si has 
been performed and the increase in melting point is explained in terms of the 
RDF peak height [15]. For Cu, at a temperature of 1440K, the intermediate 
and long distance peaks of Cu-Cu RDF are seen to slowly merge out indicating 
that melting has occured (Fig. [T]). The jump in diffusion coefficient for Cu is 
shown in Fig. [2] wherein the melting point 1340K is indicated using arrows. For 
all the doped CuTi alloys, the melting points are obtained in a similar manner. 

2.3 Embedded Atom Method Potential 

We have accounted for the interaction between atoms using the Embedded Atom 
Method (EAM) potential which incorporates the many-atom interactions ne- 
glected within the pair potential scheme. The EAM provides a semi-empirical 
potential for computing the total energy of a metallic system [T6J. In this 
scheme, we view the energy of the metal as obtained by embedding an atom 
into the local electron density provided by the remaining atoms, in addition 
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Figure 2: Jump in the diffusion coefficient (D) of Cu at melting point. 



to the pair interaction. It is especially useful for systems with large unit cells 
and is appropriate for metals with empty or filled d-bands |17j . Thus the total 
internal energy, which consists of an embedding term plus the usual two body 
interactions, is expressed as [TH] 



(1) 

(2) 



where pt is the local electron density at i atom due to all the other atoms, 
4>{rij) is the pair potential and is the distance between i th and j th atoms. 
The embedding energy, F(pi), which is a function of the local-electron density, 
is approximated by superposing the atomic electron density functions, /(r^). 
The two-body potential and electron density functions, as given by Johnson, 
are respectively written as 



(r) 



Ae- a ^-V Be-^-V 



l + (J__„)m i+(i__A)« 



(3) 
(4) 



The pair term is a Morse-type potential which is composed of a short-range 
repulsive exponential and a long-range attractive exponential. The parameters 
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in Eqs. (3)-(4) are given in Table [TJ The embedding function, F, is given by 

F(p) = F e (l-M-)")(-)' 7 (5) 

Pe Pe 

where F e is the embedding function corresponding to equilibrium electron den- 
sity p e - A smoother embedding function that matches physical properties better 
than Eq. (5) is given by [19] 

3 

F (P) = Yl F ^(P/Pn ~ 1)*) P < Pn = 0.85p e , (6) 

i=0 
3 

F (P) = E F *0>/Pe - ^ Pn < P < PO = 1.15p e> (7) 

i=0 

f( P ) = Foa-M-ni-)", pq<p (8) 

Pe Pe 

The parameters in Eqs. (6)-(8) are given in Tabled] 

For an alloy model using EAM, an embedding function F(p) and an atomic 
electron-density function / (r) are to be specified for each atomic species, while 
two-body potentials are needed for each possible pairs of atomic species. Since 
the electron density at any location is taken as a linear superposition of atomic 
electron densities, and the embedding energy is assumed to be independent 
of the source of the electron density, these two functions can be taken from 
mono atomic models directly. For a binary alloy with a- and b-type atoms, 4> aa 
and (j> bb are given by the individual mono atomic models, and <p ab and <j) ba are 
assumed to be equal. A direct and simple possible form for the alloy potential 
can assume either the geometric or arithmetic average of two-body potentials 
as has been shown by Foiles et al. [3D]. However, a novel two-body alloy model, 
which satisfies the invariance properties of electron density in a fee mono atomic 
metal, has been constructed by Johnson [21j . In this work, we adopt Johnson's 
alloy model using EAM which satisfies the invariance condition: 

where <fi ab is the alloy potential, f is the density function and the superscripts 
aa and bb stand for mono atomic a and b, respectively. 



2.4 Initial CuTi alloy structure 

Structures of CuTi alloys having different Ti concentrations have been generated 
by random or selective substitution of Ti in perfect fee Cu supercells. The origin 
of the cartesian coordinate system is taken to be the centre of the supercell. 
The perfect fee Cu supercell is generated using the program genlat.f in utility 
of DL_POLY. The Cu atoms are placed one after the other starting from the 
(-,-,-) octant towards the (+,+,+) octant. 
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Table 1: Parameters of electron density function and two-body potential for Cu 



and Ti. 





Cu 


Ti 


re(A) 


2.556 


2.930 


/ e (eV/A) 


1.554 


1.860 


a 


7.670 


8.780 


P 


4.09 


4.68 


A 


0.328 


0.328 


B 


0.469 


0.469 


K 


0.431 


0.431 


X 


0.863 


0.863 


m 


20 


20 


n 


20 


20 



Table 2: Parameters of embedding function for Cu and Ti. 
Cu Ti 



V 


0.921 


0.560 


Pe(eV/A) 


22.150 


25.600 


F„ (eV) 


-2.176 


-3.200 


F nl (eV) 


-0.140 


-0.200 


F n2 (eV) 


0.286 


0.680 


F n3 (cV) 


-1.751 


-2.320 


F (eV) 


-2.190 


-3.220 


Fi(eV) 


0.000 


0.000 


F 2 (cV) 


0.703 


0.610 


F 3 {eV) 


0.684 


-0.750 
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For random doping of Cu atoms with Ti, random numbers (depending on 
the concentration) lying within the supercells are generated for x, y and z co- 
ordinates. The Cu atom whose coordinates are nearest to the random numbers 
are replaced by Ti atoms. 

Microstructures of Ti having various sizes are also generated within the Cu 
supercells. Single microstructures of different concentrations are generated by 
replacing all the Cu atoms with Ti atoms inside spheres of varying radii at the 
centre of the microstructures i.e., (0,0,0). For substituting 8 microstructures, 
all Cu atoms lying within spheres of different radii centered at (-D,-D,-D), (- 
D,-D,+D), (-D,+D,-D), (+D,-D,-D), (-D,+D,+D), (+D,-D,+D), (+D,+D,-D) 
and (+D,+D,+D) with D=10.845A are replaced by Ti atoms. Similarly, 9 
microstructures are obtained by replacing the Cu atoms within the sphere at 
the origin in addition to the spheres used for 8 microstructures. 

Selective doping is done in a variety of ways: For a concentration of 5 % Ti, 
first Cu atom among every 20 atoms is replaced by Ti till all the 5324 atoms 
are covered. The selectively doped CuTi alloy generated in this way is called 
selective 5 % atoml. Similarly for a concentration of 10 % Ti, first Cu atom 
among every 10 atoms is replaced by Ti. This alloy is named as selective 10 % 
atoml. A type named atom2 is generated by replacing 2 atoms at a time. So, 
for a concentration of 5 % Ti, first two Cu atoms among every 40 atoms are 
replaced by Ti. Thus, selective 33.33 % atoml is generated by replacing first 
among every 3 Cu atoms with Ti and selective 33.33 % atom6 by replacing first 
six among every 18 Cu atoms with Ti. 

The natural CuTi alloy structures, viz. CuTi and CuTi2 are also generated 
using the program genlat.f. The number of atoms of Cu and Ti and their 
positions within the unit cells are obtained from ICSD database |22j . 

3 Results and Discussions 
3.1 Random doping 

The melting temperature for Cu is reported as 1200K in an earlier study em- 
ploying DL-POLY with Sutton Chen potential [14]. However, we get 1340 K 
which agrees better with the experimental result of 1356K. For random doping, 
melting point decreases linearly as Ti concentration increases (inset of Fig. [4]). 
Since the atomic radius of Ti (2 A) is higher than that of Cu (1.57 A), cell 
volume increases when a Ti atom replaces a Cu atom. The increase in lattice 
parameter of the CuTi alloy is found to obey the empirical Vegard's law, which, 
for a given temperature, is a linear relation between lattice constant and con- 
centration of the constituent elements [23] (see Fig. [3]). As a result, the average 
distance between Cu atoms increases and the Cu-Cu bond strength decreases. 
This is clearly reflected in Fig. |4]which shows the OK RDF for the Cu-Cu bond. 
It is observed that the heights of the RDF peaks decrease, full widths at half 
maxima (FWHM) increase and the peak positions shift to the right. As men- 
tioned earlier, decrease in RDF peak height is directly linked to the reduction 
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Figure 3: Linear variation of the lattice parameter as a function of the atomic 
weight percent of the dopant Ti in Cu. 

in number of nearest neighbours. Also, increase in cell volume leads to shifting 
of RDF peak to the right. Thus for random doping, the average Cu-Cu bond 
strength and hence melting point decreases linearly as dopant concentration 
increases. 

3.2 Microstructural doping 

For most of the practical cases, doping of clusters of several atoms is more prob- 
able compared to atom by atom substitution. Therefore, micro-structures of 
several Ti atoms are created within the Cu lattice . For a single microstructure, 
the Cu-Cu bond strength decreases linearly with increase in Ti concentration 
as shown in FigJSJ This is also reflected in the linear decrease of melting point 
(inset of Fig. [5J . Similarly, for 8 microstructure doping, the RDF peak height 
decreases linearly (Fig. [6]) and a gradual decrease in melting point is observed 
(see inset of Fig. |6]). Exactly similar trends are observed for 9 microstructure 
doping. 

3.3 Selective doping 

The systematic variation in melting point on changing the Ti concentration as 
observed for random and microstructural doping is no longer seen in case of 
selective dopant substitution. Both dopant concentration and the site of substi- 
tution are responsible for determining the melting point of the alloy. As already 
pointed out, Ti atoms can be introduced into Cu lattice in various ways to 
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Figure 4: RDFs for Cu-Cu bond in the case of random doping of Cu with Ti. 
Inset shows linear variation in melting point (M.P.) as a function of the number 
percent of the dopant Ti in Cu. 




Figure 5: RDFs for Cu-Cu bond in case of single microstructure doping of 
Cu with Ti. Inset shows linear variation in melting point as a function of the 
number percent of the dopant Ti in Cu. 
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Figure 6: RDFs for Cu-Cu bond in case of 8 microstructure doping of Cu with 
Ti. Inset shows linear variation in melting point as a function of the number 
percent of the dopant Ti in Cu. 



generate different atomic arrangements having the same dopant concentration. 
When a Ti atom replaces a single Cu atom in a unit (atoml), melting point de- 
creases with Ti concentration up to 20 % as in the case of random substitution 
(inset of Fig. [7]). Then the melting point increases for 25% Ti. As shown in 
Fig. Ul upto 20 % Ti concentration, the height of the first RDF peak decreases, 
become broader and position shifts to the right showing a loss in symmetry of 
the structure. But for 25% of Ti concentration, the first RDF peak becomes 
narrower and attains its maximum value. Its position does not shift further 
indicating stronger Cu-Cu bonding for the structure as compared to 20 % Ti. 
Therefore, high melting point observed in case of 25 % atoml case can be un- 
derstood in terms of the height, width and position of the RDF peak. In a 
similar manner, the low melting point found in case of 33.33 % atoml case can 
be understood by the loss in Cu-Cu bond strength which is clearly evident from 
the corresponding short and broad RDF peak. Finally, peak height increases 
and becomes narrower for 50 % Ti which results in increase in melting point. 
Similar kind of non-monotonous behaviour of melting point with dopant con- 
centration can be seen for selective atom2 to atom6 cases. Only the late rise in 
melting point can be seen for different dopant concentration depending on the 
initial atomic arrangement. 

The correlation between the melting point and characteristics of the RDF 
peak established here is important for the following reason. For an alloy it is 
necessary to quantify the interaction between the host atoms as well as between 
the host and dopant atoms. Especially, at higher temperatures when atoms 
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Figure 7: RDFs for Cu-Cu bond in case of selective doping of Cu with Ti 
(atoml). Inset shows variation in melting points as a function of the number 
percent of the dopant Ti in Cu. 

may be displaced from their equilibrium positions, characteristics of the RDF 
peaks can predict the symmetry of the atomic configuration and the thermal 
behaviour of the alloy. 

The role of site selective doping on thermal stability can be unequivocally 
established if melting point is shown to change for different arrangement of 
dopant atoms but having same concentration. For that purpose, six different 
configurations of 25% doped CuTi alloy having different arrangements of Ti 
atoms in Cu lattice are constructed. Variation of the melting point for different 
Ti arrangements is depicted in the inset of Fig. [5] The first, second and third 
Cu-Cu RDF peaks for four arrangements, namely, atoml, atom2, atom4 and 
atom6 are shown in Fig. [8] The melting point is seen to decrease upto atom4 
and then increases for atom5 and atom6. As shown in Fig. |8l up to atom4, 
the height of the first RDF peak decreases showing a loss in symmetry of the 
structure. However, for atom6, the RDF peak height increases and shifts to the 
left showing a more compact structure thus explaining the increase in melting 
point. 

Finally, we plot the melting points against dopant concentration for different 
ways of doping, namely, random, microstructural and selective doping (Fig. [§]). 
For random, single, 8 and 9 microstructures, linear variation of melting point on 
Ti concentration is observed. However, in case of microstructural doping, the 
melting points decrease faster compared to random doping. Selective doping 
(atoml), on the other hand, does not show the linear variation and the melting 
point depends on the initial structure irrespective of the concentration. 
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Figure 8: First three peaks of the Cu-Cu RDF for selective substitutional doping 
of Cu with 25% Ti doping. Inset shows variation in melting points for different 
Ti arrangements having 25% Ti concentration. 
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Figure 9: Melting points obtained for different types of substitutional doping of 
Cu with Ti. 
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Figure 10: First Cu-Cu RDF peaks for natural phases of CuTi alloy. Inset shows 
variation in melting points. 

3.4 Natural CuTi alloys 

The link established so far, between melting points and characteristics of the 
RDF peaks of an alloy, can be employed to understand the observed anomaly 
in the melting curve of naturally existing phases of CuTi alloys. Melting point 
as extracted from the phase diagram of CuTi alloy show non-monotonous de- 
pendence on Ti concentration. In the inset of Fig. [101 the calculated melting 
points of Cu and two of its naturally occurring alloys viz. CuTi and CuTi2 
are plotted and the corresponding RDF peaks are shown in Fig. [TO] In line 
with the experimental phase diagram, melting point Vs. Ti concentration first 
decreases and then again increases. For example, substantial reduction in the 
melting point of naturally occuring CuTi alloy (1210 K) is seen compared to 
pure Cu (1340 K). Melting point of CuTi2 increases to 1360K for 66.6 % doping 
(see inset of Fig. [T0| . The late rise in melting point with the increase of Ti con- 
centration arises due to higher ordering between the atoms in CuTi2 (I4/mmm) 
than CuTi (P4/nmmS) . This is reflected on the Cu-Cu bond strength which can 
be visualized by the reduced RDF peak height of CuTi compared to pure Cu. 
As expected, RDF peak height of CuTi2 increases which explains its observed 
higher melting point. 

4 Conclusions 

In summary, the role of site-selective substitution of Ti in Cu on the melting 
point of CuTi alloy has been investigated. Super cells of CuTi alloy having dif- 
ferent arrangement of Ti atoms are constructed. Results obtained by replacing 
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Cu atoms by Ti randomly, selectively and in the form of clusters are analyzed. 
We have established that, in addition to the concentration, the arrangement of 
dopant atoms in the host lattice plays a pivotal role in determining the melting 
point. A direct link between the melting point and characteristics of the RDF 
peaks of the alloy has been established. This facilitates to explain the varia- 
tion in thermal stability in terms of the bond strength between host as well 
as host and dopant atoms in the alloy. The proposition has been validated by 
explaining the anomaly in the melting curve seen in naturally occurring phases 
of CuTi alloys having different crystal structures. The present study can be 
extended to other alloys of its kind and is useful for predicting doping strategies 
for fabrication of the alloy. 
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